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A quantum theory for mesoscopic electric circuits in accord with the discreteness of electric 
charges is proposed. On the basis of the theory, Schrodinger equation for the quantum LC-design 
and L-design is solved exactly. The uncertainty relation for electric charge and current is obtained 
and a minimum uncertainty state is solved. By introducing a gauge field, a formula for persistent 
current arising from magnetic flux is obtained from a new point of view. 
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Along with the dramatic achievement in nanotcchnology, such as molecular beam epitaxy, atomic scale fabrication 
or advanced lithography, mesoscopic physics and nanoelectronics are undergoing a rapid development @^]. It has 
been a strong and definite trend in the miniaturization of integrated circuits and components towards atomic-scale 
dimensions S for the electronic device community. When transport dimension reaches a characteristic dimension, 
namely, the charge carrier inelastic coherence length, one must address not only quantum mechanical property but 
also the discreteness of electron charge. Thus a correct quantum theory is indispensable for the device physics in 
integrated circuits of nanoelectronics. Since the classical equation of motion for an electric circuit of LC-design is 
] just the same as that for a harmonic oscillator whereas the 'coordinate' has the meaning of electric charge, the 
quantization of the circuit was carried out Q in the same way as that of a harmonic oscillator. This only results 
energy quantization. In fact, a different kind of fluctuation in mesoscopic system, which inherently has nothing to do 
with energy quantization and interference of wave functions, is due to the quantization of electronic charge. Recently 
we studied the quantization of electric circuit of LC-design under consideration of the discreteness of electric charge 
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I. INTRODUCTION 



In present paper we extend the main idea of our previous letter Q and present a quantum mechanical theory 
for electric circuits based on the fact that electronic charge takes discrete values. In section |l], A finite-difference 



. Schrodinger equation for the mesoscopic electric circuit is obtained. In section [11]], the Schrodinger equation for a 



mesoscopic circuit of LC-design is turned to Mathieu equation in 'p-representation' and solved exactly. The average 



value of electric current for ground state is calculated. In section [V, the uncertainty relation for charge and current 
' is discussed, a minimum uncertainty state, which recovers the usual Gaussian wave packet in the limit of vanishing 
, of discreteness, is solved. In section |v|, the Schrodinger equation for a quantum L-design both in the presence of an 
adiabatic power source and in the absence of source are solved exactly. A gauge field is introduced and a formula 



for persistent current which is a periodic function of the magnetic flux is obtained. It provides a formulation of the 
persistent current in mesoscopic ring from a new point of view. Finally, some discussions and conclusions are made 
in section |Vl]. 



II. QUANTIZATION OF ELECTRIC CIRCUIT IN ACCORD WITH THE DISCRETENESS OF 

ELECTRIC CHARGE 

We recall that for a classical non-dissipative electric circuit of LC-design in the presence of a source e(t), the 

d 2 q 1 1 

equation of motion, as a consequence of Kirchoff's law, reads H -pzq s(t) = 0, where q(t) stands for electric 

charge; L for inductance and C for the capacity of the circuit. This equation of motion can be formulated in terms of 
Hamiltonian mechanics, namely 
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. _ dH dH 

^ dp ' P dq 

with H(t) = — p 2 H — — q 2 + e(t)q. Here the variable q stands for the electric charge instead of the conventional 

'coordinate', while its conjugation variable p(t) = L— represents (apart from a factor L ) the electric current instead 

of the conventional 'momentum'. Analogous to the forced harmonic oscillator, the electric circuit was quantized by 
many authors Q, where the electric charge was treated as a continuous variable. As a matter of fact, the electronic 
charge is discrete and it must play an important role in the theory for mesoscopic circuits. Taking account of the 
discreteness of electric charge, we must reconsider the quantization of a mesoscopic circuit. According to the standard 
quantization principle, one associate with each of the two observable quantities q and p a linear Hermitian operator, 

namely q and p. The Hamiltonian, also an observable quantity, corresponds to a Hermitian operator H = -zyp 2 + V(q) 
which is a function of the operator p and q. The commutation relation for the conjugation variables are 

[q,p]=ih (2.1) 

Up to now, the discreteness of electronic charge is not taken into account. Regarding to the discreteness, we must 
impose that the eigenvalues of the self-adjoint operator q take discrete values j5), i.e. 

q\q>=nq e \q> (2.2) 

where n € Z (set of integers) and q e — 1.602 x 10 -19 coulomb, the elementary electric charge. Obviously, any 
eigenstate of q can be specified by an integer. This allows us to introduce a minimum 'shift operator' Q :— e lqe ^^ h , 
which is shown to satisfy the following commutation relations § 

[<i, Q] = -QeQ 

[q 1 Q + ]=q e Q + 

Q+Q = QQ+ = 1. (2.3) 
These relations can determine the structure of the whole Fock space. For q\n >= nq e \n >, the algebraic relations 



(2.3) enable us to derive the followings 

Q+\n > = e iQ "+ 1 |n+ 1 >, 
Q\n >= e- ian \n- 1 >, (2.4) 

where q„'s are undetermined phases. Obviously Q + and Q are ladder operators respectively for charge increasing 
and decreasing in the diagonal representation of charge operator. The Fock space for our present algebra differs 
from the well known Fock space for the Heisenberg-Weyl algebra, because the spectrum of the former is isomorphic 
to the set of integers Z but that of the later is isomorphic to the set of non-negative integers Z + + {0}. Since 
{\n > 1 7i 6 Z} constitute a Hilbert space, we have the completeness J^ntz \ n >< n \ = We also have the 
orthogonality < n\m >= S nm due to the self-adjointness of q. As a result, the inner product in charge representation 
takes as 

< <f>\ip >= V" < <f>\n >< n\ip >= V" <f>* (n)i[>(n) (2.5) 
ne z ne z 

One can now study the eigenstates and eigenvalues of the operator p. Obviously, if p\p >= p\p > then f(p)\p >= 
f(p)\p > for any analytical function /. Supposing \p >= J2 neZ c n (p)\n >, and using Q\p >= e lq " p l h \p > we can find 
that Cn+i/c-n = <2xp(iq e p/h + ia n +i), which yields the following solution 

\p>=J2 K ne inq " p/h \n> (2.6) 

n£ z 

where n n — e 1 ^-^ 1 ^, k_„ = e *S^=o a ^ 3 for n > 0. Obviously \p + h(2n/q e ) >= \p >, the eigenvalues of the 
operator p is a periodic parameter. Topologically, the parameter space of the spectrum is isotopic to the S 1 . 
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Since the spectrum of charge is discrete and the inner product in charge representation is a sum instead of the usual 
integral, one may define a right and left discrete derivative operators V ge and V gc by 



V,./(») = 
V ?e /(n) = 



f(n + l)-f(n) 
/(n)-/(n-l) 



(2.7) 



They can be understood as the inverse of a discrete definite integral, which is in accord with the inner product (|2.E|) , 



/ f{x)dx := ^2 q e f{nq e ) 



__ f QF(xf) - F{xi) ifV q .F = f, 
\ F(x f )-Q+F{ Xi ) ifV qe F = f, 

Clearly, it recovers the conventional diffe rent ial-integral calculus as long as the minimum interval q e goes to zero. The 
discrete derivative operators defined by (2.7) can be expressed explicitly by the minimum shift operators 



V 



(Q - l)/«e 

(1 - Q + )/q. 



(2. 



It is easy to check H that V 
operator 



V 9c . Then we can write down two important self-adjoint operators: 'momentum' 
K ,_ = s h 



^=2i(V,.+V 9a ) 



2iq, 



(Q-Q + ) 



and free Hamiltonian operator 



■(v,.-v g .),= 



2<?i 



■(Q + Q+-2) 



(2.9) 



(2.10) 



where we call them respectively as momentum and free Hamiltonian operators because they are really those when 
q e — * 0. Now we have finished the quantization of mesoscopic electric circuits and obtained the following finite- 
difference Schrodinger equation, 



2q e L 



(V„.-V,.) + V($) 



\ip >= E\tp > . 



(2.11) 



III. THE QUANTUM LC-DESIGN 



As an application of our quantization strategy of the mesoscopic circuit, we discuss a mesoscopic LC-design in this 
section. We on ly consider the adiabatic approximation so that e(t) is consider as a constant e. Then the Schrodinger 
equation (2.11) for a LC-design is written as 



2q e L 



(V, 



V 9 J + ^ 



eq 



\iP >= E\ip > 



(3.1) 



We consider a representation in which the operator p is diagonal and called it as p-representation. We must 
addr ess that the p is the conjugation of the charge variable q within the meaning of usual canonical commutator 
( [2.l|) , and it is the 'current' operator only if the charge is treated as a continuous variable. However, the operator 

P associated with the physical quantity, electric current(apart from a factor — ), differs from the operator p as long 



as the discreteness of charge is taken into account. Clearly, P will become the usual p when q e goes to zero. The 
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orthogonality of eigenstates of p is an immediate consequence of (2.6) and the orthogonality of the charge eigenstates, 
i.e. < p\p' >= — v 5(p — p' + n( — )h). The completeness is also verified 



h{ t) dp 
2tt J jv_\ h 



q, 



p >< p\ = \n >< n\ = 1. 



- n (^r) ■■ nez 

The transformation of wave functions between charge representation and p-representation is given by 

< n\ip >= ( — -) / dp < p\ip > e m b 
2nh J-n(f-) 



(3.2) 



(3.3) 



Using (2.6), we can obtain the following relations 



< p \ q \p > 



4irk 
2nh 3 d 2 



< P'|V ge - V, e \p > = -^r (cos(|p) - l) S(p - p') 

2 S(p-p') 



q e dp 2 



(3.4) 



In the 'p-representation', the finite-difference Schrodinger equation ( |3.l| ) becomes a differential equation for ip(p) :=< 
P\ip > 



2C dp 2 q 2 L 



4>(p) = E4>(p). 



(3.5) 



which is the well known Mathieu equation [|7|Jg] . This equation was ever appeared in [Bl on the discussion of Pade 
approximates. In deriving to (3.5), we have adopted e = for simplicity. Actually, the linear term in ( |3.l| ) can be 
moved by a translation i n th e 'coordinate' (charge) space. Apart from a re-definition of q and a shift of the energy 
E, the same equation as (3.5) would be derived. 

In terms of the conventional notations , the wave functions in p-representation can be solved as follows 

4>+(p) = ce/(| - f) 



TT q e 



(3.6) 



2~h C 

where the superscripts '+' and '-' specify the even and odd parity solutions respectively; I = 0, 1, 2, • • •; £ = (^r) 2 — ; 

qi L 

ce{z,£) and se(z,£) are periodic Mathieu functions. In this case, there exist infinitely many eigenvalues {a;} and 
{bi+i} which are not identically equal to zero. Then the energy spectrum is expressed in terms of the eigenvalues a;, 
bi of Mathieu equation 



E, 



i+i 



Si 

8C~ 



h 2 
qlL 



(3.7) 



As an exercise, one may calculated the fluctuation of electric current for the ground state. It is known that the 
explicit results of eigenvalues and eigenfunctions of Mathieu equation are complicated. They are related to continue 
fractions and trigonometric series respectively. For the concrete values of the Plank constant and the elementary 
electric charge, the WKB method is valid. From the series solution of Mathieu equation for ground state, we obtained 
the fluctuation of electric current P (apart from a factor 1/L) for ground state, 



< P 2 > A= - I — 

\ j .--ground I 



_ ^ ( %2c 

~2 Up 



(3. 



This result is valid for the case C/L << {q 2 /T%y 
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IV. UNCERTAINTY RELATION AND THE MINIMUM UNCERTAINTY STATE 



In order to fairly understand the main conclusions in this section, we begin with a brief view of the derivation of 
the uncertainty relation in standard quantum mechanics. If A and B are two Hermitian (self-adjoint) operators which 
do not commute, the physical quantities A and B cannot both be sharply defined simultaneously. The variances of 
A and B are defined as (AA) 2 =< (A— < A >) 2 > and (AB) 2 =< (B— < B >) 2 >. Their positive square roots, 
AA and AB are called the uncertainties in A and B. In terms of the properties of self-adjoint operators and the 
knowledge of Schwarz inequality, one can prove that 



(aA) 2 (ab) 2 



>|<^({i,B}- 



< A >< B >) > 



\<\[A,B}>\ 2 



(4.1) 



where { , } denotes the anti-commutator, and the equality sign holds if and only if B\ip >oc A\ip >. In deriving (4.1), 
the fact that the expectation value of a Hermitian (or an ti-Hermitian) operator is a real number (or purely imaginary 
number) has been used. As a direct consequence of (4.1), the uncertainty relation is conventionally written as 



{aA) 2 {ab)'- 



>\<^ B}>\ 2 



(4.2) 



Clearly, the equality sign in (4.2) holds if and only if both the equality sign in (4.1) holds and the first term of the 
right hand side in (4.1) vanishes. These conditions imply that 



(B— < B >)\ip > 



A = 



A(i- < A >)\ip > 
<#i, B]\jy> 

2(Ai) 2 



(4.3) 



Now we direct to our main purpose. After some calculations, we obtain the following commutation relations for 
the charge q, the current P and the free Hamiltonian Hq 



[Ho, P] = 0, [d , q] = ihP, [q, P] = iTl(l + f§ff ) 

n 



(4.4) 



where the operators P and Hq have been defined respectively by (2.9) and (2.10). The term ^\Ho in the third 

h 



equation of (4.4) occurs due to the discreteness of electric charge. Now we are ready to write out the uncertainty 
relation for electric charge and electric current, namely 



Aq-AP> J(i + || < JT >). 



(4.5) 



This uncertainty relation recovers the usual Heisenberg uncertainty relatio n if q e goes to zero, i.e. the case that the 
discreteness of electric charge vanishes. Moreover, the uncertainty relation (4.5) has shown us some new sense beyond 
the knowledge of traditional Heisenberg uncertainty relation. 

It is of interest to study the particular state \4> >, for which (JT^) becomes an equality. This is the state in 
which the product of the uncertainties in electric charge and current is as small as the noncommutivity allows: 



Aq-AP=-(l- 



— |- < Ho >)• Such a minimum uncertainty state must obey the condition (4^) for A = P and B = q 



{q-<q >)|V >= 



ih(l + -% < H >) 



2(AP) 2 



-{P- < P >)\iP > 



(4.6) 



Using (^6) and (2^), one can find that 



<P q P> = — -p-<5(p-p) 
q e % op 

< P '\P\p> =--sm(^f)S(p-p') 

q e q e n 



(4.7) 
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Then (4.6) becomes the following differential equation in p-representation 

„2 



2(AP) 2 Ue R ' 



(4.8) 



This differential equation is solved by a plane wave with modulated amplitude: 



ip(p) = N exp 



<H > , %2 



2(AP) 2 



h " f1 e \ , ^ 6 ^ \ i<q>p 
— cos(— p) + < P > p ) - 

9e 



(4.9) 



where N is the normalization constant. ( [4.9] ) is obviously a deformation of the usual Gaussian wave packet and 
recovers the Gaussian wave-packet if the discreteness vanishes. 



V. QUANTUM L-DESIGN, GAUGE FIELD AND PERSISTENT CURRENTS 

In this section, we will solve the Schrodinger equation for a L-design in the presence of an adiabatic source and in 
the absence of source. Introducing a gauge field and gauge transformation, we derive a formula for persistent current 
in a pure L-design, i.e. a mesoscopic metal ring. 



A. The L-design in the presence of an adiabatic source 



The Schrodinger equation for a L-design in the presence of an adiabatic source reads 



2q e L 



(V,.-V,.)+ea 



\ip >= E\ip > 



(5.1) 



In order that the quantization of a mesoscopic circuit be valid, the size of the circuit must be restricted. While 
the voltage source can come from an infinite reservoir to keep the chemical potential constant. We consider present 
problem in charge representation and expand the eigenstate of (5.1) in terms of orthonormal set of charge eigenstates, 
namely, \ip >= X) n =-oo Un \ n - > - Substituting it into (5.1), we obtain the following recursion relations 



/ h 2 k 2 . 

2( — r + nq e e - E)ui - - — -(uj_i + Uj+i) = 
q e L 2q e L 



(5.2) 



The knowle dge of recursion formula of Bessel functions, z(J„ + i(z) + J v -\{z) ) — 2vJ v (z) enables us to write down 
a solution of (5.2) 



J, 



nq e s+zo — E 



(*b) 



(5.3) 



where zq = . In terms of (5.3) the eigenstates of Schrodinger equation (p. if) are write out 

q e L 



\lpE >= ^2 Jnq a e+z -E(zo)\n > 
n— — oc 

which is the solution of eigenstates for quantum L-design in presence of a adiabatic source. 



(5.4) 
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B. Pure L-design 



h 2 

Now we consider a pure L-design, Hl = — — — p(V ge — V ge ). The Hamiltonian operator of the pure L-design Hl, 

2q e L 

t he c urrent operator P and the operator p commute each other, so they can have simultaneous eigenstates. Actually, 
(2.6) is the simultaneous eigenstate of those operators, i.e. 

P\p > = — sm( — ) \p >, 

Qe ft* 



ffib> =| z (l-cos(^)) |p>, (5.5) 



This result tells us that the magnitude of electric current in a mesoscopic electric circuit of pure L design are bounded 

taking values between — ( ) and ( ). It also indicates that the maximum quantum noise in a pure L-design 

q e L q e L 

( a mesoscopic ring is an example) takes finite value if the elementary charge q e should not be considered as the 
infinitesimal ( particularly for the mesoscopic circuit). It is also worthwhile to notice that both the current and 

energy of a pure L-design become null when p = as long as q e is not zero. Clearly, the lowest energy states 

'/■ 

correspond to p = n — for any integer n. Thus the energy spectrum is infinitely degenerated. 
q e 



C. Gauge field and persistent current 

In previous discussion, we have used the terminology 'p-representation' and solved the eigenstates of p. Now let us 
to find out what the eigenvalues of p mean. If introducing a operator G :— e~ l ^i , we can find that G\p >= \p — (3 > 
and G + \p >= \p + (3 >. Considering a unitary transformation to the eigenstates of Schrodinger operator given by 

|V> >^ W >= G\1> >, 



we find that the Schrodinger equation (2.11) is not covariant. This requests us to introduce a gauge field and to define 
a reasonable covariant discrete derivative. By making the following definitions: 



GD^G- 1 = D' qe , 



we can verify that they are covariant under a gauge transformation. The gauge transformations are expressed as 

v-i _ ( 

GD^G- 1 =T)' qe , (5.7) 
as long as the gauge field <p transforms in such a way 

^ <t>' = <l> - 13. 

From either the transformation law or the dimension of the field </>, we may realize that <f> plays the role of the magnetic 
flux threading the circuit. 

In terms of those covariant discrete derivatives (5.6), one can write down the Schrodinger equation in the presence 
of the gauge field (magnetic flux). Here we write out the Schrodinger equation for a pure L-design in the presence of 
magnetic flux, 

-^L^s -DM>=E\tP>. (5.8) 
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£(p,0) = ^sin 2 (fl( P -0)) (5.9) 



Because its eigenstates can be simultaneous eigenstates of p. ( p.8[ ) is solved by the same eigenstate \p > in (2.6). The 
energy spectrum is easily calculated as 

2h . 2 ( ' Qe 

H " "' 

which has oscillatory property with respect to <j) or p. Differing from the usual classical pure L-design, the energy of 
a mesoscopic quantum pure L-design can not be large than 2fijq\. Clearly, the lowest energy states are those states 

that p = <f> + n( — ). Thus the eigenvalues of the electric current ( i.e. — P ) of ground state are calculated 
qe L 

J(^) = A flin (^). (5.10) 

Obviously, the electric current on a mesoscopic circuit of pure L-design is not null in the presence of a magnetic 
flux except <fi = n( — ). Clearly, this is a pure quantum characteristic. (5.1C) exhibits that the persistent current in 

Qe 

a mesoscopic L-design is an observable quanti ty pe riodically depending on the flux cj). Because a mesoscopic metal 



ring is a natural pure L-design, the formula (5.10) is valid for persistent current on a single mesoscopic ring [ p_0| 
Differing from the conventional formulation of the persistent current on the basis of quantum dynamics for electrons, 
our formulation presented a method from a new point of view. Formally, the I(<j)) we obtained here is a sine function 

h 



with periodicity of 4>q = — , But either the model that the electrons move freely in an ideal ring 1 1 1 , or the model that 

^ e • in 

the electrons have hard-core interactions between them [n2| can only give the sawtooth-type periodicity. Obviously, 

the sawtooth- type function is only the limit case for q e /h — * 0. 

Certainly the experiment Jl3| should be considered as the case of persistent current in a LC-design because the 
junction of semiconductors will contribute a capacitance to the 'circuit'. 

VI. CONCLUSIONS AND DISCUSSIONS 

In the above, we studied the quantization of mesoscopic electric circuit. Differing from the literature in which it 
is simply treated as the quantization of a harmonic oscillator, we addressed the importance of the discreteness of 
electric charge. Taking the discreteness into account, we proposed a quantum theory for mesoscopic electric circuit 
and give a finite-difference Schrodinger equation for mesoscopic electric circuit. As the Schrodinger equation for 
LC-design in p-representation becomes the well known Mathieu equation, it is exactly solved. We obtain the wave 
functions in terms of Mathieu functions and the energy spectrum in terms of the eigenvalues of Mathieu equation. 
The discussion on uncertainty relation for the charge and current shed some new light on the knowledge of transitional 
Heisenberg uncertainty relation. The discreteness of electric charge increased the uncertainty which is related to the 
expectation value of the 'free' Hamiltonian. The minimum uncertainty state we obtained is a deformation of the 
standard Gaussian wave packet. As further applications of our theory, the eigenstates of L-design in the presence and 
in the absence of source were solved respectively. Introducing a gauge field and gauge transformation, we successfully 
obtained a formula for the persistent current on the mesoscopic pure L-design in the presence of the magnetic flux. 
As the mesoscopic metal ring is a natural pure L-design, the formula is certainly valid for the persistent current on 
mesoscopic rings. In our formula, the mass of electrons, the carriers for electric current, is not involved. This is 
worthwhile to check by some experiment. Our present theory is believed to explain the Coulomb blockade on which 
the research is in progress. 

In addition, all the results in present paper will recover the standard knowledge if one takes the continuous limit 



q e — > 0, e.g. ( 5.10 ) becomes 4> = LI in the limit of q e — > 0, the well known formula in electromagnetism. So the 
whole theory and their results are believed to be consistent and reasonable. One may noticed that we used the charge 
representation and the so called p-representation. Because of the discreteness of electric charge, p is no long er a 



current operator, but should be understood as the usual Dirac conjugation of the charge operator satisfying (2.1) 



only. The operator P which is associated with physical observable, electric current, obeys the commutation relation 



(4.4). Clearly, the current operator P is not a Dirac conjugation of the charge operator q. So we need a new definition 



about such conjugation defined by (4.4). 
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